Reconstruction of f(R) models with Scale-invariant Power Spectrum 
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Following our previous work in [JCAP 1206, 041 (2012) [l[], in this paper, we continue our study 
of reconstructing f(R) modified gravity models that can be connected to a single scalar field in 
general relativity via conformal transformation, which lead to scale-invariant power spectrum in the 
early universe. With f(R) modified gravity, one does not need to introduce extra scalar, the nature 
of which are to be explained. Different from general nonminimal coupling theory, the behavior of 
the f(R) theory has been fixed by its counterpart in Einstein frame, and thus have one to one 
correspondence. Numerical plots of the functional form of f(R) as well as the evolution of R in 
terms of cosmic time t are also presented. 



I. INTRODUCTION 



For theories of the early universe, the right amount of 
perturbations must be generated so as to conform with 
our observations such as cosmic miscrowave background 
(CMB) Q and large scale structures (LSS) Q. One of the 
well-known observed features is, the power spectrum of 
these perturbations, which comes from the 2-point cor- 
relation function, has to be (nearly) scale-invariant 
which will put on non-trivial constraints on theoretical 
model building. Although it is well-known that a single 
scalar field, which drives the universe into de-Sitter like 
expansion (inflation [4-8], while the scalar is called in- 
flaton), or nonrelativistic matter-like contraction |qL Hoj 
could easily generate perturbations to meet the require- 
ment, the nature of the scalar is still unclear. 

Scale-invariant power spectrum may also arise when 
one modify Einstein's gravity at early times. In some 
cases, the modified gravity theories could be connected 
with unmodified general relativity (GR) plus a scalar 
through conformal transformations with the latter 
being viewed as the counterpart in Einstein frame of the 
former. Due to the equivalence between the two frames 
(Jordan and Einstein) , the perturbation generated by the 
couple of counterparts are exactly the same. Thanks to 
the connection, one can thus reconstruct models of modi- 
fied gravity from the known evolution of GR plus a scalar 
models, which can lead to inflation or matter-contraction 
scenarios. Recently we proposed a way of reconstructing 
the models with a scalar nonminimally coupled to gravity 
which could give rise to scale- invariant power spectrum 
In this paper, we will consider another case of mod- 
ified gravity, namely f(R) theories. Actually as we will 
see later, f(R) theories could be one specific but non- 
trivial form of nonminimal coupling. In f(R) theories, 
there is no need to introduce the unknown scalar, and 
the universe is driven totally by its gravitational struc- 
ture. f(R) theories has been used widely as alternatives 
of inflation, dark matter, dark energy and so on. See [l2[ 



for comprehensive reviews. 

The reconstruction of f(R) gravity has been pursued 
by many authors, see [13]. In their approaches, most 
of them reconstruct f(R) theory in Jordan frame itself, 
provided that the cosmic evolution in Jordan frame is 
given. Here we will reconstruct in a different way, namely 
from their counterpart in Einstein frame, which looks like 
a single scalar field in GR, via conformal transformation. 
This kind of reconstruction aims at connecting different 
evolutions of the universe driven by modified gravity in 
its Jordan and Einstein frames. As is shown in [l[, in 
Einstein frame there are only two cases which could give 
rise to (nearly) scale-invariant power spectrum, namely 
inflation and matter-contraction. Taking the Einstein 
frame lagrangian as: 
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Re - AdifE? - V(cp E ) 



(1) 



where here and after we set the unit such that 8irG = 
Mpf = 1, and use the metric signature (—,+,+,+). A 
simple and representative solution is the exact solution 
which is obtained assuming that its equation of state w E 
is a constant, namely: 



a E (t E ) ~ (±^)3(i+™ E > , H E (t E ) 



3(1 + w E )t E 



<PE(t E ) = 



2H±Mt E ) T „.„_, , e -V5(I+=5)* 



^3(1 +w E ) 



f2) 



where M is some energy scale. In this parametrization, 
we have set " + " for positive t E meaning an expanding 
phase, while " — " for negative i_g denoting a contracting 
phase, and Vo is some constant factor. In Inflation case, 
we have w E = — 1 + 2e E /3 with the slow- roll parameter 
\e E \ = I - (dH E /dt E )/H 2 E \ < 1, then Eq. © can be 
written as: 

aE(t E ) 



t E 'E , H E (t E ) 



£ E t E 



VE{t E ) = J— HMt E ) , V(v E ) = V Q e-^* E (3) 

V e E 
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while in matter-contraction case, one has w E — 0, and 
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thus Eq. becomes: 



Perturbations 



a E {t E ) ~ (~t E ) 3 , ffs^B) 



3t f 



<p E (t E ) = _ln(-A/t £ ) , V{ip E ) = V e-^ E .(4) 

In this short paper, we will mainly focus on the Jor- 
dan frame of the modified gravity theories in order to find 
which form can be conformally connected to the above 
two cases, while more complete study for the case of vary- 
ing w E (or e E ) will be left for the future. 

The remaining sections are organized as following: in 
Sec. II, we review the main results for the general non- 
minimal coupling theories that was obtained in our pre- 
vious paper, in Sec. Ill, we focus on f(R) theories. Nu- 
merical plots of the functional form of f(R) as well as the 
evolution of R in terms of cosmic time t are presented. 
Furthermore, we also discussed about the relation of the 
evolutions of various cosmological variables between the 
two frames for an arbitrary constant e E . In Sec. IV we 
conclude our paper. 



II. REVIEW OF RECONSTRUCTION OF 
NONMINIMAL COUPLING THEORY 

A. Background 



The equation of motion of the perturbation generated 
by the action ([5]) can be written down as: 



aj\/2Qiz 



)«R = , 



(9) 



where u-ji — a,jy/2Q-ji7Z, and 1Z is the conformal- 
invariant curvature perturbation. The variable Q-n is 
defined as: 



2F 



(2 + 8 F ) 2 



(10) 



where 8f = F/(HjF). The prime denotes derivative 
with respect to the conformal time rj — J aj 1 (tj)dtj. 
By solving Eq. ©, it is straightforward to see that 
scale-invariant spectrum can be obtained in two ways, 
one is a,jy / 2Qn ~ \rj* — rj\ where the time- varying 
mode becomes decaying while the constant mode domi- 
nates the perturbation, which is inflation, and the other 
is ajy / 2Qiz ~ |?y* — rj\ 2 where the time- varying mode is 
the growing mode and thus dominates over the constant 
one, which is matter-contraction. In fact, from the rela- 
tion ([5]) one can express Qn as: 



(11) 



First of all, we will briefly review the main results ob- 
tained in [![. The action of the nonminimal coupling 
theory we are considering is: 



and since we have assumed constant w E and e E , the con- 
dition of getting scale-invariant power-spectrum can be 
written as aj\fF ~ — ?7 1 1 or a,j Vf~ 



s 



NMC 



1 



(5) 

where F((j>) and Z((j>) can be arbitrary functions of the 
field 4> in the Jordan frame, and £/(</>) is the potential. 
The equation of motion of is: 



3Hj(f> 



2Z 



(H + 2H 2 ) 



z 



0, (6) 



where subscript "0" indicates d/d(j) and dot denotes 
derivative with respect to cosmic time tj in the Jordan 
frame, and the Friedmann Equation is: 



QHjF + 6HjF = \z$ + U 



(7) 



Following the conformal transformation of metrics in Jor- 
dan and Einstein frame, gffi = Q 2 g$ , where fl 2 = 2F, 
the relations of some basic variables between the two 
frames are summarized as follows: 



dt i 



Qdt 



J , a E 



H 



"(I 



<p E = 



Z 



O 2 



-d<j> , V(cp E ) 



2Hjtt' 

u(4>) 



n 4 



(8) 



C. Reconstruction of nonminimal coupling theory 
in Jordan Frame 

We can reconstruct the universe evolution once we as- 
sume the evolution of in terms of tj. In our previous 
paper [l|, we assumed that f2(tj) = f2o[(±tj)/(±tj)]", 
then from the relation (El) we have: 



t E 



n+i \±t*j) 



for uj 5^ — 1 



(12) 



n t*j\n(±tj) for ui = -l 



where the " + " sign in " ± " means tj > 0, and in the 
Jordan frame the universe is expanding, while the " — " 
sign means t j < 0, and in the Jordan frame the universe 
is contracting. Here we define tj = tj/tpi where tpi is 
the Planck time. Substituting it into Eqs. ([3]) and (|4]) 
respectively, one can get the evolution of variables such 
as aj, Hj and wj in terms of tj as (for u> ^ — 1 only): 



aj(tj) ~ (±tjY 



Hj 



l + (l-e E )u 
e E tj 



Wj 



2 

-1 + - 



c E 



31 + (l-e E )w 



(13) 
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where \e E \ <C 1 for the case corresponding to inflation, 
while e E = 3(1 + w E )/2 = 3/2 for the case correspond- 
ing to matter-contraction, respectively. Moreover, from 
relation ([5]) one can also find the evolution of field vari- 
ables, and thus determine the form of functions F((f)), 
Z(<fr) and U(4>) in the lagrangian. In fact, taking the 
ansatz of Z(<fi) — Zo4> 2z and U((f>) = Uo4> q , and with the 
help of Eqs. ^ and ([7]), we found the relation: 



F(4>) =F ^ Z+Z ,q 



2(z + 1)(1- ~), (14) 

and the equation of state wj can be given by: 

_ 2(z + l)(5e g -6)-g(2e g -3) 
WJ - 3[2( Z +l)(2-e E )-q) ' (15j 

From above we can see that, once the functional form 
of F(<p), Z{4>) and U(<j>) in action (0 is given by (Tl4l) . 
one could obtain scale-invariant power spectrum. Rather 
than being fixed to be inflation or matter-contraction 
only, the evolution of the universe in the Jordan frame 
has more freedom. This is because in the Jordan frame, 
the nonminimal coupling action ([5]) has more degrees of 
freedom than that in the Einstein frame and is more de- 
pendent on the form of the action. However, as we will 
see below, it is not the case in f(R) theory. In f{R) 
theory, there will be less degree of freedom than non- 
minimal coupling theory and the form of f(R) will be 
more fixed. Following similar steps, we will find the ap- 
propriate f{R) theory, which can correspond to inflation 
or matter-contraction scenarios in its Einstein frame and 
thus, give rise to scale-invariant power spectrum. 



III. RECONSTRUCTION OF f(R) MODIFIED 
GRAVITY THEORY 



respectively. 

Same as nonminimal coupling theory, f(R) theories 
with action (|16|) can also be connected with (TTJ) as its 
counterpart in the Einstein frame, via the conformal 
transformation gff = £l 2 g\?J with il 2 — 2F. To see 
this, one can rewrite the action (|16|) in the form of scalar- 
tensor theory, namely as: 



-g F(R)R~U{R) 



(19) 



where the potential U(R) can be identified as F(R)R — 
f(R). The relations of the basic variables between the 
two frames are summarized as follows: 



dt E = VLdtj , a E = ttaj , H E = ^-—(1 , 

Si ZHjit 



ip E = v61nQ , V(ip E ) 



u{R) 
si 4 



(20) 



From the transformed action ([T9|) we can see that, the 
f(R) action is actually the specific form of the general 
nonminimal coupling action ([5]) with Z((f>) = 0, as long 
as we identify F(cf>) with F{R), and U(<fi) with U(R), 
which is easy provided that the inverse function of F((f>) 
exists. Moreover, since Z{4>) as well as the kinetic term 
of (0 vanishes, there are less degrees of freedom in f(R) 
than in nonminimal coupling theories, and the confor- 
mal factor SI, which determines the cosmic evolution 
in Jordan frame, can be totally fixed by the field <p E . 
Therefore, when there is one kind of evolution in Ein- 
stein frame, there is only one kind of evolution in Jordan 
frame. This gives less possibilities for f(R) theories to 
get scale- invariant power spectrum than those for non- 
minimal coupling theories. 



A. Background 

Now we turn on to study the reconstruction of f(R) 
modified gravity theories. The action of f(R) modified 
gravity theory is: 



(16) 



where f(R) can be arbitrary function of the Ricci scalar 
R. Varying the action (fTo]) with respect to the metric 
9lii> we can get the equation of motion: 



1 



g^OF + FR^- = 



(17) 



where we defined the function F(R) = df/dR. The left 
part of the above equation can also be viewed as the "ef- 
fective" stress energy tensor of f(R) modified grav- 
ity, which satisfies the continuity equation, V'E^ = 0. 
Moreover, the "0 — 0" and "0 — i" components of Eq. (fl7|) 
are just Friedmann equations, which are 

3H 2 F =^(f + 3F + 3HF) , — 2HF = F — HF , (18) 



B. Perturbations 

One can also check from the perturbation theory of 
f(R) that what conditions should be met when one re- 
quires a scale-invariant power spectrum. Working in the 
Arnowitt-Deser-Misner (ADM) formalism [lj], one can 
obtain the perturbed action of f{R) up to the second 
order as: 



5 (2) = / dr]d 3 xa 2 jQn \n' 2 - (dK) 



(21) 



where 1Z is the conformal-invariant curvature perturba- 
tion, and 



6F6 2 F 
(2 + S F ) 2 



(22) 



with Sp = F/{HjF) and the prime denotes derivative 
with respect to the conformal time 77. Varying (|2ip with 
respect to 1Z, one can straightforwardly write down the 
equation of motion for the perturbation as: 



(k 



2 (OJVTO) 



ajy/2Q n 



-)u n = , 



(23) 



4 



through the redefined variables u-ji — ajy/2Q-jzR-- 

From the above analysis, we can directly conclude that 
scale-invariant spectrum can be obtained in two ways, 
| — vl^ 1 which corresponds to in- 
\f]* — i]\ 2 which corresponds to 



namely aj^/2Q-jz 



flation, or ajy/2Q-ji 
matter-contraction. Moreover, from the relation (|20|) 
one can express Q-r as Q-ji ~ Fee, the same as that 
in nonminimal coupling theories. Here we can see again 
that f(R) theories are nothing but specific case of non- 
minimal coupling theories. In our case where constant 
we and ee have been assumed, the condition of get- 
ting scale-invariant power-spectrum can be written as 
ajVF — |j7» - rj\~ l or ajVF ~ I??* - t]\ 2 . 



C. Reconstruction of f(R) modified gravity theory 
in Jordan Frame 

First of all, from relations (f2T))) as well as the evolution 
of pe^e) in the Einstein frame ([2]), we can obtain the 
evolution of the conformal factor f2 in terms of ts, which 
is 



where t*, 



n 



M~ 



ft_E_ 

\t* 



\e E \ < 1 



(24) 



Since the universe in Einstein frame 

Since 



is expanding, we set ts and t* E to be positive [IE 
dtj = fl~ 1 (tE)dtE, one could easily get tj as: 



tj = - 
or equivalently, 

t_E_ _ f -tj 

t% ~ v-n 



3 ( t E 

tvh: 



(25) 



(26) 



Note that since \ee\ <C 1, tj and t*j < 0. Then we have: 

n(t J ) = [^)^~ 1 . (27) 

With Eqs. ©, (!| and (|27j) in hand, we can obtain 
the evolution of aj, Hj and wj in the Jordan frame, in 
terms of tj. The results are: 



frame) is less than zero, and aj(tj) will be increasing 
as tj increases. This indicates that it is an expanding 
universe, driven by f(R) modified gravity theory, which 
is equivalent to the so-called "Super-inflation" [l5[ (or 
phantom- inflation [l6j]) scenario in GR when transformed 
to the Einstein frame. One can also look into the equa- 
tion of state wj of the universe, which is very much close 
to —1 up to order of slow-roll parameter, which means 
that the universe in the Jordan frame is also near de Sit- 
ter, so different from the general nonminimal coupling 
theory, inflation in the Einstein frame can only refer to 
inflation in the Jordan frame in f(R) modified gravity 
theory. 

The Ricci scalar R, which is defined as R — 6(H + 
2H 2 ), can be expressed as: 



R(tj) = 6 



(2- V3ii)(3-e E ) 



(29) 



Finally, with Eqs. ([27]), (HU), as well as the relation fi 2 
2F, we can obtain the form of F(R) as: 



and 



f(R) 



R™ f = 6 (2 



/3ej)(3 - (e) 



(30) 



F(R)dR 



--R 



in J 



\ R "ifJ 



(31) 



We can see that when ce is small during inflation, the 
function of f(R) is almost proportional to R 2 up to slow- 
roll parameter. Therefore, this model coincides with the 
well-known Starobinsky's model [|| of which f(R) ~ R + 
aR 2 in the very early time, when R is very large. In the 
late time when ee is large, it goes near the standard GR. 

The plot of R(tj) and f(R), which are reconstructed 
from inflation in its Einstein frame, are presented in Figs. 
□ andU 



The Ricci scalar R(t_J) reconstructed from inflation. 

R(U) 



aj{tj) 
Hj(tj) 



Wj 



(^) 

V3 



V3- 



Te + 2\/3e_E — 3v3 



3(V3 - y/IE) 



(28) 



From this result we can see that, since \ee\ "C 1 as we 
considered, the index of aj in terms of tj (namely 1/ej, if 
we define ej to be the slow-roll parameter in the Jordan 



FIG. 1: The behavior of R(tj) w.r.t. tj, where we choose 
M = 0.1 and hence t* E = 10. In this case, R > 0, and 
is increasing w.r.t. tj, showing a "super/phantom-inflation" 
behavior. 

Following the same procedure, we can do the recon- 
struction of f(R) from matter contraction, just replacing 



5 



f(R) reconstructed from inflation. 




FIG. 2: The behavior of f(R) w.r.t. R, where we choose M = 
0.1 and hence t* E — 10. We can see that f(R) monotonically 
increases with R, and in the limit of large R, it approaches to 
the squared power- law f{R) ~ R 2 - 



tE by — ts, and e# by the value 3/2. Note that here i# 
and t* E are negative, tj and £l{tj) will become 



-tE 



,0 = y(3 



E ■ 



and 



2 + 3V2" 



(32) 



(33) 



where tj and t*j still smaller than 0. The scale factor aj, 
the Hubble parameter Hj and the equation of state wj 
will be given by: 



aj(tj) 



-tj 

~t*r 



4—V2 



, Hj(tj) 



4-V2 1 + V2 

-nT' WJ = ^- 

(34) 

From this result we can see that, since the index of 
aj in terms of tj is larger than zero, so aj(tj) will be 
decreasing as tj increases, indicating that there is also 
an contracting universe driven by f(R) modified grav- 
ity theory when we require it be equivalent to matter- 
contraction scenario in GR when transformed to the Ein- 
stein frame. The Hubble parameter Hj{tj) is smaller 
than zero because of the negative tj, and the equation of 
state wj of the universe is about the value of 0.8, which 
is even larger. 

The Ricci scalar R in this case is: 



R(tj) = 



6(8 - 9^2) 



which gives the form of F(R) as: 



2 + 3s/2 



R: 



MC 



_ 6(8 - 9V2) 



49i* 



(35) 



(36) 



and 



f(R) 



F(R)dR 



5 + 3\/2 



R 



MC 



(r* ic ) 



(5 + 3%/2) 



DA 

>Rq 



(37) 



The plot of R(tj) and f{R), which are reconstructed from 
matter-contraction in its Einstein frame, are presented in 
Figs. [3]and |U 

The Ricci scalar R(t_J) reconstructed from matter-contraction. 

R(U> 



FIG. 3: The behavior of R(tj) w.r.t. tj, where we choose 
M = 0.1 and hence t% = 10. In this case, R < 0, and is 
decreasing w.r.t. tj. 



f(R) reconstructed from matter-contraction. 




FIG. 4: The behavior of f(R) w.r.t. R, where we choose 
M — 0.1 and hence t* E = 10. We can see that f(R) is also 
less than 0, and increases with R since both R and f(R) is 
decreasing w.r.t. tj. 

One can check our results with the conditions for gen- 
erating scale-invariant power spectrum for consistency. 
For the case of inflation, from Eq. (|2"51) we can write 
down the relation of conformal time r\ and tj as: 



V = I a J 1 (tj)dt J 



(-tj)- 



while 



%/3 



(38) 



(39) 



where we note that 5f is a constant. Thus we could easily 
find that 



l/{e B -l) 



(40) 



when I eg I -C 1. The case of matter-contraction is similar. 
From Eq. (|3~4")l one has: 



/ s 3+ VI 



and 



2(3+V2) 



(41) 
(42) 
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is 




a E 


(= l/V3e B ) 


(~ - w)Jt B ) 


/ ( 1 \ // , 1/ \\ 

(= (<" - - 1/es)) 


(~*y ) 


horizon 
problem 


t E > o 


es > 3 


expanding 


1/eg < oj < 1 


tj > 


ej < 


contracting 


V 


1 < EE < 3 


w < 1/e.E 


ej > 1 


expanding 


1/3 < e E < 1 


w < 1 


< ej < 1 


n 


< e B < 1/3 


1 < U) < l/t E 


tj < 


ej < 


t E < 


€E > 3 


contracting 


1/tE < OJ < 1 


tj < 


ej < 


expanding 


n 


1< e B < 3 


U < 1/tE 


ej > 1 


contracting 


1/3 < e E < 1 


U < 1 


< ej < 1 


y 


< £E < 1/3 


1 < w < 1/es 


tj > 


ej < 



TABLE I: The relations between variables in the Jordan and Einstein frames where we generalize e_g to be an arbitrary positive 
constant value, ts can be chosen as either positive or negative, presenting parametrization of an expanding or a contracting 
universe. For eg > 1 in expanding phase or e_g < 1 in contracting phase, we have horizon problem, while in the other two cases 
we don't. Due to the fact that oj — l/v3es, the region of tj > 0/ < can be divided by the line of cj = l(eB = 1/3), the region 
of ej > 0/ < is divided by both ui = 1 and uj — 1/e.E (ee = 1). In the ej > region, the region of ej > 1/ < 1 is divided by 
the line of tj — 1/e.E (e.E = 3). Whether the universe contracts or expands in the Jordan frame is decided by whether e^ > 3 or 
not. Finally, when there is no horizon problem in the Einstein frame, there will be no horizon problem in the Jordan frame, 
and vice versa. Similar summary but only for GR can be found in, e.g. p^| . 



which gives 

ajVQ^~ v 2 . (43) 

Before ending this section, let's also remark the rela- 
tion between the general evolutions of the universe driven 



by f(R) modified gravity theory in the two frames with 
an arbitrary constant eg, though without showing the 
detailed calculations. The relation between variables in 
the two frames is summarized in TABLE [TJ Here we write 
£l{tE) in a general form of il — (±te/ ± t* E ) u . 



Note that since f(R) theory can only be equivalent to 
canonical field via conformal transformation, we don't 
have 6e < case. 



IV. DISCUSSIONS AND CONCLUSION 

In this paper we studied the reconstruction and cosmic 
evolutions of f(R) modified gravity models, which can be 
transformed as inflation or matter-contraction scenarios 
in their Einstein frame. The equivalence of the Jordan 
and Einstein frames guarantee that the perturbations 
generated by f(R) models follows the same evolution, 
namely can give rise to scale-invariant power spectrum re- 
quired by observations, however their background evolu- 
tion might be different. In our previous work [1| we have 
shown that there can be more than one kind of evolution 
in the case of general nonminimal coupling theories, but 
for the f(R) case, there's no such degeneracy and the cor- 
respondence between the two frames must be one to one. 
We find that in f(R) modified gravity theory, inflation 
in the Einstein frame can only refer to (phantom-like) 
inflation in the Jordan frame, while matter-contraction 
in the Einstein frame can only refer to contraction with 
a larger equation of state in the Jordan frame. We anal- 
ysed the general conditions for f(R) theory of getting 
scale-invariant power spectrum, and obtained the evolu- 
tion of the universe in the Jordan frame as well as the 
functional form of f(R). Numerical plot of R w.r.t. tj 
and f(R) w.r.t. R are also presented. 



As a side remark, we would like to mention that due 
to the equivalence of the two frames, the Big-Bang cos- 
mological problems (horizon, flatness, etc.) will also do 
no harm to the reconstructed f(R) models. To see this, 
one can look into the efolding number j\f defined as [l8j 



TV = In 



O-iHi 
CLiH; 



(44) 



which can be directly related to these problems. Usually 
these problems can be avoided as long as we require that 
j\f > 70 during inflation. From the relation ([20)) we can 
see that the conformal Hubble parameter, H = aH, is 
not conformal invariant, but since in our case ^ is a 
constant, Af is a conformal invariant variable. Therefore, 
provided that inflation lasts for enough efolding number 
in one frame, one need not worry about whether it does 
in the other frame. We'd also like to refer the readers to 
for more detailed arguments. 
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